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^ ' Abstract 

B: 

I Earlier Arnold, Zeldovich, Ruzmaikin and Sokoloff [JETP (1982)] have computed the 

^ ■ eigenvalue of a uniform stretching torus transformation which result on the first Riemann 

(N : 

Tj- , metric solution of the dynamo action problem. Recently some other attempts to obtain 

O I Riemann metrics representing dynamo action through conformal maps have been undertaken 

^ ■ [Phys. Plasmas 14 (2007)]. Earlier, Gilbert [Proc. Roy. Soc. London A(2003)] 

O ■ has investigated a more realistic dynamo map solution than the one presented by Arnold 

.!_H I et al by producing a shearing of the Arnold's cat map, by eigenvalue problem of a dynamo 

X ■ 

' operator. In this paper, the eigenvalue of the Riemann twisted torus dynamo flow metric is 

obtained as the ratio between the poloidal and toroidal components of the flow. This result 
is obtained from the Euler equation. In the twisted torus, the eigenvalue of the Riemann 
metric is m± = , which is quite close to the value obtained by Arnold. In the case the 
viscosity Reynolds number Re oo, the torus flow is unstable as one approaches the torus 
axis. In Arnold's dynamo metric the eigenvalues are x± = ^^2^ which are very close to the 
above value. Eigenvalues determine the growth rates of the velocity ratio between poloidal 
and toroidal components of the flow. The curved flow in torus follow previous work by Chen 
et al [Phys Fluids (2006)]. The a-effect dynamo is shown to be a second-order effect in the 
torus curvature and velocity flow. Loop dynamo flows and maps are also discussed. PACS 
numbers: 02.40.Hw:differential geometries. 91. 25. Cw-dynamo theories. 
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I Introduction 



The investigation of rotating torus flows and their perturbation [1] have been recently linked to 
the dynamo experiments with liquid sodium and other liquid metals [2, 3]. The incompressible 
flows, where most of the anti-dynamo theorems are applied [4], can be used as rotational flows 
in Perm torus [2] and Riga experiment [3]. Therefore a good understanding of the behaviour of 
these rotating torus flows is of utmost importance in planning these experiments. Alternative 
dynamo experiments using plasma flows have been recently obtained by Wang and his group 
[5], which has developed the first flowing magnetic plasma (FMP) experiment, called P — 24 
in order to detect dynamo action. In this paper a slight variation of the usual Riemann-flat 
(Riemann curvature zero) torus named a twisted torus [6] filled with a rotating incompressible 
rotating flow. The incompressibility of the flow is well suited for liquid metals. In this twisted 
torus, as shown here, the Riemann curvature tensor does not vanishes and for thick tubes 
the torus is also stretched, a fundamental property for existence of dynamo action [7]. This 
kind of twisted geometry is very usual in plasma devices such as heliotrons, stellarators and 
in the astrophysical and solar flux tubes [8]. Previously Mikhailovskii [9], have made use of 
non-diagonal Riemann metrics to describe such a plasma devices. Here one follows Ricca [8] 
approach and use a diagonal simpler twisted magnetic flux tube Riemann metric. The first 
example of a chaotic fast dynamo solution was found by Arnold et al [10] by making use of 
a compressed and stretched Riemannian metric of the dynamo flow. The domain of the flow 
is given by a compact Riemannian manifold represented as the product of a torus and the 
closed interval [0, 1] of < z < 1 [11]. This results in the Arnold Riemann metric 

ds'' = e-^'dp'' + e^'dq" + dz'' (I.l) 

which represents the stretching and contraction in distinct Euclidean directions p and q , re- 
spectively induced by the eigenvalues x± = also corresponding to magnetic eigenvectors. 
The stationary dynamo flow considered by Arnold was given by the simple uniformly stretching 
flow V = (0,0,^) in {p,q,z) coordinates. More recent attempts to build a dynamo action by 
making use of compact Riemannian geometry includes the case of the fast dynamo of Chiconne 
and Latushkin [12], and the conformally stretched fast dynamo by Garcia de Andrade [13]. 
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The rotating torus flow equations are 

Vb 

dt\i = -2VtXM + v/\m (1.2) 

Po 

where A := V^. The incompressible flows are described by the solenoidal vector field as 

V.v = (1.3) 

Equation (1.2) is the Navier- Stokes equation for the rotating flow inside the torus, while v is 
the viscosity constant. Here 

dt = dt + v.V (1.4) 

Note that the eigenvalue problem in this case is not easy to define since even in the absence 
of pressure gradients the eigenvalue would be zero and if one tries to use the ABC flow for 

example 

O = Av (1.5) 

the first term on (1.2) vanishes but yet the eigenvalue problem is not helped much. In this paper 
however, one presents an eigenvalue like problem that solves the differential Euler equation 
(1.2). The scalar proportionality between the poloidal and toroidal flows modulo the eigenvalue 
m shows that this eigenvalue is quite similar to the torus eigenvalue obtained by Arnold. 
Actually from the mathematical point of view these transformations define an automorphism 
on the torus of the type ^ : v — > v. This paper is organised as follows. In section II the 
Euler equation free of viscosity in the background of a twisted torus flow. In this case the 
radial pressure is unstable closer to the internal torus axis, as happens in some torus plasma 
instability associated to coronal mass ejection, recently investigated by Toerok and Kliem [14]. 
The a — dynamo effect is also discussed in section III. Section IV addresses the loop filaments 
dynamo torus twist maps, while discussions and conclusions appear in section V. 
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II Stretch-twist dynamo maps 

Following previous work by Childress and Bayly [13], Gilbert [15] developed applied their 
stretch-fold-shear (SFS) dynamo nnap in Arnold's cat map, by performing matrix transforma- 
tions to his the expression of cat dynamo Jacobian matrix [4] 

Mcat= (^^ 1^ (n.6) 

corresponding to the map x — > Mx. these cat maps with shear, built by Gilbert result on the 
class of matrices in the form: 

Meat-shear = ^ J (11-7) 

where K is an integer related to the eigenvalues. In this paper one shows that the torus twist 
map [6] 

which is clearly distinct from the chaotic dynamo map matrix (11.6), can be associated to the 
Ricca's Riemann twisted magnetic flux tube metric 

= dr^ + r^dOn^ + K^{r, s)ds^ (II.9) 

under the transformation 

■"t-t^li ^(!)) <"-^°* 

which is obtained from the transformation x — Or and y — s, which from the twist transfor- 
mation 9{s) := Or ~ J T{s)ds, results in matrix (11.10), as long as one considers the helical 
case where the Frenet torsion r(s) = tq = constant. In the case of thin twisted magnetic 
flux tube, K — 1 and the equation simplifies to 

Mthin-tube= ^"-^^^ 

which is still closer with the twist torus map. If one normalizes the torsion tq = —1 the last 
matrix reduces exactly to twist torus map matrix (11.8). Of course the presence of stretching 
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factor K in the Riemann metric can be considering as a constant stretch which reduces the 
eigenvalue nnatrix to 

Mt„ist-(; ";°) (11.12) 
Let b is the transpose of an initially constant toroidal magnetic field given by 

btwist = (0 1) (11-13) 

Within these motivations, in the next section, one shall follow a less mathematical path to 
help experimental dynamo physicists not very familiar heavy framework of dynamo theory, and 
compute the solution of the Euler equation with rotation in the case of a homothetic motion 
between the poloidal and toroidal flow velocities. It is shown that the self-induction magnetic 
equations eigenvalues imply stretching or contraction of the magnetic flows and growth or 
decay of the toroidal field with respect of the poloidal field, in either case dynamo action is 
present. 



Ill Unsteady twisted torus dynamo flows 

Let us impose the constraint d^v = on the Euler equation 

dtV = -2Qxv (III. 14) 

Solutions of these equations yields naturally the ABC flows above. The Euler equation can 
now be written in the background of the above twisted flux tube Riemann metric, where 
9{s) := Or — J T{s)ds and tq is the constant radius of the constant cross-section flux tube, 
and K{s) = (1 — rK,{s)cos6{s)). If the tube is thin factor K(s) pa 1 and the gradient operator 
is given by 

V ^tK-^ds + ee-de + Brdr (III.15) 

r 

where the Riemannian line element is in general given by 

dp = Qijdx'dx^ (III. 16) 

where (i,j = 1,2,3) and e R^. To express Euler equations in this metric background one 
must use the dynamical relations from vector analysis and the theory of curves in the Frenet 
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frame (t, n, b) are 

t' = Kn (III.17) 
n' = -«;t + rb (111.18) 
b' = -rn (III.19) 
The dynamical evolution equations in terms of time yields 

t = [n'h - urn] (III.20) 

li = KTt (111.21) 

b = -n't (III.22) 

along with the flow derivative 

i = dtt+ {vV)t (III.23) 

With these mathematical tools in hands now one is able to write down the Euler equations if 
one expresses the velocity flow in the form totally confined inside the torus as 

\i^Vs{r)t + ve{r,s)e0 (III.24) 

where Vs and vq are the toroidal and poloidal flow components. The Euler equations yield 

drP 2 

— = VsKq^ - veujo + - [vs- voIkqVs (III.25) 
Po r 

one notes that this first equation is non-linear on velocity flow and Kq is the Frenet curvature. 
Besides one has considered that the pressure p — p{r). The remaining equations are 

^Vs + ^v'e + v"0^jV0 (III.26) 

where the dash represents the radial derivatives. The remaining equation is 

-iv0-Vs) + -v's + v"s = -fVs (III.27) 
r r 

where the compact Riemannian operator [17] 
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and the self-induction or dynamo equation can be written as 

C^B = 7B (111.29) 

To obtain these equations a strong hypothesis simplification of 9 « 1 of small angle ampli- 
tudes is undertaken. Here one also assumes that the poloidal rotation loq inside the torus is 
constant and vq = cuor. These unsteady flow equations on the twisted torus are easily solved 
by assuming the eigenvalue like expression 

V0 — mvs (III. 30) 

where m e R and the terms non-linear in the Frenet curvature were dropped. By taking the 
coordinate substitution r' = Inr the equations and the compact Riemannian operator can be 
simplified since 

By substituting these constraints on the equations (111.26) and (111.28) yields the equations 

2vs + mv' s + mv''^ = ^mvs (III. 32) 

where the dash represents the radial derivatives. The remaining equation is 

{m-l)vs + v's + v'\^-fVs (III.33) 

Multiplying the last equation by m and subtracting the result from the first equation (111.30) 
one obtains 

[2 - m(m - l)]vs = (III.34) 

This equation can be expressed in the eigenvalue matrix form 

/(2-m(m-l)) 0\ 
Mdyn=( Q J (ni.35) 

which yields the simple eigenvalue equation 

- m - 1 = (III.36) 

which, in turn, yields 

m± = (III.37) 
7 



which is similar to the above Arnold's eigenvalue for the uniform stretching Riemann dynamo 
metric. Note that if one chooses m+ the poloidal velocity is higher than the toroidal one while 
when m_ is chosen the toroidal velocity flow dominates in intensity over the poloidal one. The 
last constraint on the torus flow equations is the solenoidal vorticity Q 

V.n = (III.38) 

yields 

v's = -- (III.39) 

r 

which yields Vg — —Inr. Substitution of this expression into the pressure equation 

p{r) = PqIujo't — mKQ{lnrY] (III. 40) 

Note that as r — > this reduces to 

p(0) = — pomKo{lnry — > oo (III. 41) 

which shows that the pressure "explodes" at the internal torus axis. This shows that the flow 
is unstable inside the torus in the absence of dissipation. One is lucky enough to understand 
that in laboratories this is a very unlikely situation since as far as liquid metals and even laminar 
plasmas are concerned the magnetic and fluid Reynolds numbers are finite and as one shall 
see in the next section in these cases the flow is stable inside the torus, which turns the task 
of building a liquid sodium dynamo flow detector in laboratory less difficulty. 
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Let us now use the computation of the vorticity used above and given by 

1 

n = --[tcosO - hsece]KoVs (111.42) 
to compute the a-dynamo [16] effect 

a = n.v (III.43) 



which then yields 



a ^ ^[m - 1]ko'^Vs^ (III.44) 



which shows clearly that the a effect is a second-order effect in the torus curvature. Substi- 
tution of the m values, after a straightforward and simple computation one obtains 

a=-^-^^^^KoW (III.45) 

from which one notes that at the torus internal axis the a effect is quite strong enhancing the 
dynamo effect. Of course the a expression is actually random and this cannot be generally 
true. 



IV Loops in twist torus non-dynamo maps and flows 

In this section one shall use the mathematical formalism of the Frenet frame of the first section 
to build the loop twist torus dynamo filamentary flows using the eigenvalue problem. Let us 

consider the matrix of the magnetic field in the Frenet frame in the form 

B=(5,0,0) (IV.46) 

and the self-induction equation 

dtB = B.Vv + jyAB (IV.47) 
By taking the compact Riemannian dynamo operator, without stretching terms as 

£^ = ryV^ - 7 (IV.48) 

One shall note that due to the absence of the stretch term one shall have to constraint direction 
of the rotational flow vorticity in order to have a non-dynamo, since as one shall note from 
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the computations the Frenet torsion shall vanish and therefore the flow is a planar flow and 
since incompressibility is assumed, from the Zeldovich anti-dynamo theorem dynamo action 
shall not be possible. By making use of the Riemannian line element along the filament 

ds^o = Ko^ds'^ (IV.49) 

This metric can be obtained from the torus or flux tube metric above, by taking the filament 
or very thin tube approximation of r ~ 0. By Taking into account that now the gradient 
operator is written as 

V = tKo-^ds (IV.50) 

Substitution of (IV.46) in the form B — Bt and the gradient operator (IV.49) into equation 
(IV.47) leads to the matrix form 

fl + ^A \ 
Mfii - 7I = -^o~'7 ^ (IV.51) 

Here A :— Kqk'k, B = C = ^Vq. Where the 3x3 matrix missing term M33 equals 

-tKq but due to the computations this term vanishes due to the vanishing of torsion. Here 
(1 0\ 

1=1 j is the unit matrix. Taking the determinant of the expression (IV.51) one 
obtains the eigenvalue equation as 

i:>ei[Mfi| - 7I] = (IV.52) 

which implies the algebraic equation 

^{BAC) + C+{^yBA^Q (IV.53) 
7 7 

which results in 

^{BAC)+C+{^fBA^O (IV.54) 
7 7 

which indicates that either the dynamo is slow, since 7 may vanishes when the diffusion 
r] vanishes, or the Frenet curvature is constant. In this last case the dynamo could be fast, 
however this is an unphysical solution since in turbulent flow for example, a constant filamentary 
curvature should be unstable and variations in curvatute along the filament would produce 
again a slow or marginal dynamo. Actually, the vanishing of filament torsion implies that the 
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flow is planar and since incompressibility condition was used throughout the computations 
according to Zeldovich anti-dynanno theorem there is no dynamo action at all. Actually this 
was expected due to the absence of the stretching terms. Some maps as the Baker's maps 
are non-dynamo maps [18, 19] such as this one one just considered here. 
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V Conclusions 



In this paper the importance of investigation the twisted torus dynamo unsteady flows, in 
non-dissipative case for the building of the torus dynamo experiments with liquid metals was 
discussed. In particular it is shown that dynamo flow solutions can be found with special 
eigenvalue problem. More recent chaotic fast dynamos were obtained by Lau-Finn result 
which has a growth rate of 7 = 0.077 for Rtm = 00. They also obtain for Rem = 1000 
a 7 = 0.076. Growth rates as high as 3 have been found in Perm liquid sodium dynamo 
torus, as has been computed by Dobler et al [2]. The ideas developed here followed the path 
of realistic dynamo maps obtained from the Gilbert's shearing generalization of Arnold's cat 
map. Topology of chaos [20] is a nice subject which deals with torsion and knott theory 
, writhe numbers to useful in the further investigation of the topics discussed here. In the 
filamentary case discussed above actually an application of anti-dynamo theorem of Vishik's 
[21] has been done since slow dynamos in non-stretching flows are obtained. 
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